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ABSTRACT 

The light curve, energy spectra, energetics, and IPN localization of an exceedingly 
intense short duration hard spectrum burst, GRB 070201, obtained from Konus-Wind, 
INTEGRAL (SPI-ACS), and MESSENGER data are presented. The total fluence of 



x 



the burst and the peak flux are S = 2.00+q;26 x 10 ~ 5 er S cm ~ 2 and -^max = 1-611J 
10~ 3 erg cm~ 2 s _1 . The IPN error box has an area of 446 square arcminutes and covers 
the peripheral part of the M31 galaxy. Assuming that the source of the burst is indeed in 
M31 at a distance of 0.78 Mpc, the measured values of the fluence S and maximum flux 
-Pmax correspond to a total energy of Q = 1.5 x 10 45 erg, and a maximum luminosity 
L = 1.2 X 10 47 erg s _1 . These data are in good agreement with the corresponding 
characteristics of the previously observed giant flares from other soft gamma repeaters. 
The evidence for the identification of this event as a giant flare from a soft gamma 
repeater in the M31 galaxy is presented. 

Subject headings: gamma-ray bursts; soft gamma-ray repeaters 



1. INTRODUCTION 

Soft gamma repeaters (SGRs) are a special rare class of strongly magnetized neutron stars ex- 
hibiting two types of gamma-ray burst emission. Occasionally, SGRs enter an active stage and emit 
repeated short bursts with spectra which in an energy range above 10-15 keV can be approximated 
by a soft thermal bremsstrahlung-like function with kT ~ 20-30 keV. The bursting activity may 
last from several days to a year or more, followed by a long quiescent period lasting up to several 
years. 

Much more rarely, perhaps once every 30-40 years, an SGR may emit a giant flare (GF) with 
enormous intensity. The energy released in such a flare in gamma-rays is comparable to the total 
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energy emitted by the Sun over 10 4 -10 5 years or more. GFs display short (~ 0.2-0.5 s) initial 
pulses of hard gamma rays with a steep leading edge and a more extended trailing edge which as 
a rule evolves into a soft long decaying tail pulsating with the neut ron star rotation perio d. A GF 
is frequently preceded by pronounced increase in bursting activity (jFrederiks et al.ll2007al ). 



The first two SGRs were discovered and localized in March, 1979 by the Konus experiment on 
Venera 11 and 12. A giant flare from SGR 0526-66 was detected o n March 5, 1979 by numerous 
spacecraft (jMazets et al.lll979al : Evans et al.lll980l : ICline et alJll980l ). One of the first results of the 
interplanetary network (IPN) was the localization of the fi rst GF to an erro r box inside the N49 
supernova remnant in the Large Magellanic Cloud (LMC) (jCline et al.l ll980). Another important 
result was the observation and localization by the Konus e xperiment of ~17 weak soft recurrent 



burst s emitted by SGR 0526-66 in the years following the GF (jMazets et al.lll979al : lGolenetskii et al 



19841 ). Moreover, in 1979 Marc h, Konus detected an d localized three similar recurrent bursts from 
anoth er source, SGR 1900+14 dMazets et al.lfl979bh. The third S GR, 1806-20, was discovered in 
1983 dLaros et al.lll987j: ktteia et al.lll987l : Kouveliotou et al.l ll987n. and the fourth, SGR 1627-41, 

1999bl : 



in 1998 (H urlev et al 



Woods et al 



1999). It has become clear that SGRs are very rare 



astrophysical objects. The giant flare on 1979 March 5 remained a unique event for 20 y ears, until 
new giant flares were detected from SGR 1627-41 on 199 8 June 18 (jMazets et al.lll999al) and from 



SGR 1900+14 on 1998 August 27 (jHurlev et al. 



1999a 



Feroci et al 



1999; 



Mazets et al 



Final l y, a giant flare from SGR 1806-20 on December 27, 2004 was also observed (jHurlev et al 



2005 



1999b). 



Frederiks et al.ll2007al ). The extremely high intensity of the initial pulses of the giant flares 



suggested that these events could be detected in galaxies at distances up to 10-30 Mpc. 

The short hard GRB 051103 was localized by the interplanetar y network to the nearb y M81 
group of interacting galaxies which lies at a distance of 3.6 Mpc (jGolenetskii et al.ll2005l ). The 
possibility of identifying this burst as a GF from a n SGR in the M81 group, as opposed to a short 
GRB in a background galaxy, has been discussed (jFrederiks et al.ll2007bl ; iHurlev et al.ll2007al ). 



On 2007 February 1 Konus- Wind recorded a short duration hard spectrum gamma-ray burst 
whic h had a higher intensity than any previously observed event. The burst was localized by the 
IPN ( Golenetskii et al.1 120071 ). The final localization (a 446 arcmin 2 error box) covers the outer 
arms of the M31 galaxy. The analysis presented below of the time history, spectral characteristics, 
and energetics argues strongly for the interpretation that this ev ent is actually a GF from a soft 
gamma repeater in the Andromeda galaxy. In a companion paper. lAbbott et al.l (|2007l ) discuss the 
LIGO data taken at the time of this event. 



2. OBSERVATIONS AND LOCALIZATION 

The light curve of GRB 070201 (Konus trigger time To=55390.780 s UT, corresponding to an 
Earth-crossing time 55390.261 s) recorded by the Konus- Wind S2 detector in the energy range 17- 
1130 keV is shown in Fig. [U It displays a narrow pulse with a rather steep leading edge (~20 ms) 
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followed by more prolonged decay to T-To ~ 180 ms. The maximum count rate occurs in a < 2 ms 
long interval. The burst profile is not smooth: in most cases, the variations evident in the count rate 
at the top of the pulse and in its decay phase are statistically significant. A weak secondary flash 
is observed at T-To ~ 180-280 ms. Comparison of the light curves in three energy bands Gl (17— 
70 keV), G2 (70-300 keV), and G3 (300-1130 keV) (Fig. [2]) and of the hardness ratios G2/G1 and 
G3/G2 reveals a strong, rapid spectral evolution of the radiation. The most pronounced changes 
are observed for the high-energy part of the spectrum (energies > 300 keV) . The duration of the 
high-energy radiation does not exceed 80 ms. The second, weak rise in intensity at T-To ~ 180 ms 
is present only for the soft part of the spectrum at energies E^ < 300 keV. Four multichannel 
spectra were measured in the course of the burst, each with an accumulation time of 64 ms. They 
cover the energy range from 17 keV to 14 MeV, but no statistically significant emission is seen 
above 2 MeV. The boundaries of the successive accumulation intervals N=l, 2, 3, and 4 are shown 
in Fig. [2] by dashed vertical lines. Because of the small number of counts in interval 3, spectra 3 
and 4 were combined. The time-integrated spectrum of the entire burst was accumulated over the 
interval 0-256 ms. The raw count rate spectra were rebinn ed in order to have at least 20 counts 



per energy bin and then fitted using XSPEC, version 11.3 (|Arnaudl ll996). The detector response 
function was calculated specifically for the 55?9 burst incidence angle, determined from the IPN 
localization data. 

A good fit was obtained for a power-law spectrum with an exponential cutoff, dNp^/dE cx 
E~ a exp [—(2 — a)E/E p ], where E p is the peak energy in the EF{E) spectrum. We also tested the 
Band (GRRM) model but no statistically significant high energy power-law tail was found in any 
fitted spectrum. The deconvolved photon spectra for the intervals 1, 2, and 3+4, and for the time- 
integrated spectrum of the burst, are shown in Fig. [3l and the spectral parameters are summarized 
in Table [H Comparison of the E p values shows that the spectrum for interval 1 is the hardest. 
Returning to Fig. [21 note that the accumulation time of this spectrum is rather long (64 ms) 
compared to the characteristic evolution time and hence the spectrum is strongly time-averaged. 
Most of the high-energy photons are accumulated during the first ~ 20 ms. This remarkable fact 
means that the flux of photons with the hardest spectrum is emitted near the peak of the initial 
pulse. The burst fluence in the 20 keV-1.2 MeV range is 2.00] l Ig'2g x 10~ 5 er g cm -2 . The 2-ms peak 
flux measured from To+0.016 s in the same energy band is 1.61±g;|g x 10 erg cm s . All the 
uncertainties are for the 90% confidence level. 



GRB 070201 was detected by Wind (Konus) with a time resolution of up to 2 ms, by INTE- 
GRAL (SPI-ACS) with a resolution of 50 ms, by MESSENGER (MErcury Surfac e, Space ENvi- 



ronme nt, GEochemistry, and Ranging , Gamma-Ray and Neutron Spectrometer, iGoldsten et al 



(|2007l )) with a resolution of 1 s, and by the Swift BAT with 1 s resolution (outside the coded 
field of view, J. Cum mings, private communication 2007). The initial triangulation was given in 
Hurley et al.1 ( 2007bl ). The coordinates of the center and vertices of the refined 446 arcmin 2 3a 



error box are listed in Table El The center lies ~ 1° away from the center of M31. Assuming that 
the source of GRB 070201 is situated in M31 at a distance of 0.78 Mpc, the measured values of the 
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fluence and peak flux correspond to an isotropic energy output Q = 1.5 x 10 erg, and an isotropic 
peak luminosity L = 1.2 x 10 47 erg s _1 . 

The relative proximity of M31 makes it worthwhile to search for the afterglow of GRB 070201 in 
soft gamma rays, i.e., the tail of the possible giant flare. While recording a burst in the trigger mode 
with high time resolution, Konus also continues to measure the count rates in each of the detectors 
SI and S2 in four energy windows Gl, G2, G3, and Z (a charged particles channel) in background 
mode with a time resolution of 2.944 s. The initial pulse of GRB 070201 falls completely in a single 
2.944 s interval around To- We began by selecting a long series of count rate data, excluding this 
interval. Then, in order to suppress fast statistical fluctuations we formed a new series consisting 
of a sum over 94.2 s (32 x 2.944s), which can originate at the beginning of any 2.944 s interval. Two 
such series, in the energy bands Gl and G2, are shown in Figures [4^, and b; they are synchronized 
to the first 2.944 s interval following the excluded one. The position of the excluded interval is 
marked by a narrow vertical line. The discontinuity at T-To ~240-3800 s is due to the transfer of 
the accumulated data to the on-board memory. A significant increase in the count rate (4.3<r above 
the average background level) is seen only in the soft energy band Gl and only for the first large 
interval after To. Panel shows how the excess counts AiV vary if the beginning of the 94.2 s 
interval is shifted forwards or backwards from To step-by-step by 2.944 s. The dependence of AiV 
on the time shift AT implies that the gradually decreasing soft gamma-ray flux is indeed present in 
each 2.944 s interval after To but is absent prior to it. This important result confirms the detection 
of the decaying soft afterglow of GRB 070201. For an OTTB spectrum with kT ~ 30 keV, the 
94.2 s count excess AiV = 1385 ± 320 corresponds to a fluence S ~ 1 x 10~ 6 erg cm" 2 . Assuming 
that we are observing emission in the tail of a GF from M31, its energy is Qtaii ~ 7 x 10 43 erg. 



3. DISCUSSION 



The Andromeda galaxy, M31, as the closest massive galaxy, has long been regarded as one 
of the most likely candidates for searching f or observable GFs from distant extragalactic SGRs 
(jMazets et al.lll999bl ; iBisnovatvi-Koganl 11999) . It is well known that one of the domin ant f e ature s 
of the structure of M31 is a bright circular ring with a radius of 10 kpc discovered by lArpI (j 19641 ) 
in Ha observations. In agreement with numerous observations, the ring reveals the main star- 
formation region of the galaxy. In the far infrared, the first IRAS images of M31 also revealed a 
bright circular ring (jHabing et al.lll984l ). The close correspondence between the d etails of the Ha 



and f ar infrared images of the galaxy implies that the radiation of massive hot stars (jDevereux et al 



19941 ) is the common energy source of these emissions. Recently, high spatial resolution images of 
M31 have been obtained in the mid- and far infrar ed at different wavelengths with the Spitzer space 



telescope ( B armb vet al 



2006 



Gordon et al.ll2006l ) and in the far and near ultraviolet with GAL EX 



(IThilker et al.ll2005l ). These very similar images clearly show the central region, the fragmented 
spiral arms, the main circular ring with R ~ 10 kpc, and a less luminous outer ring with R ~ 14 kpc. 



Figure [5] shows the UV image of M31 (IThilker et al.ll2005l ). The IPN error box overlaps the 
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northeastern part of the galaxy. It contains segments of both the R ~ 10 kpc ring and of the 
weaker R ~ 14 kpc outer ring. The probability of a chance overlapping of the error box and 
the i mage is rather low ( ~ 10 -4 ). Twenty-two sources from an XMM-Newton X-ray survey of 



M31 (jPietsch et al. 1120051 ) fall within the error box. Among these, two sources are active galaxy 



nuclei (AGN) and five are foreground objects. SGRs in the quiescent state are typically sources of 
pulsating soft X-rays with an average 0.5-10 keV luminosity of 10 35 -10 36 er g s -1 and a spectrum 



which can be described as a combination of a blackbody and power law (IKulkarni et al 



Hurley et al 



2000 



Kouveliotou et al. 



2001 



2003 



Spectral studies of the unidentified soft X-ray sources 
in the error box which could be SGRs in M31 would clearly be important. 

We now examine additional evidence which strengthens the interpretation of the 2007 Febru- 
ary 1 event as a GF from an SGR in M31. This comes from a comparison of the energetics, time 
history, and spectral behavior of GRB 070201 with the corresponding characteristics of the previ- 
ously observed giant flares on 1979 March 5 (SGR 0525-66), 1998 June 18 (SGR 1627-41), 1998 
August 27 (SGR 1900+14), and 2004 December 27 (SGR 1806-20). Here we are mainly interested 
in the characteristics of the GF initial pulses, but unfortunately, it is just this information which 
is the hardest to obtain reliably. When a sensitive gamma-ray detector records an initial pulse 
with an immense intensity, it is strongly overloaded, which makes it difficult or even impossible to 
obtain information about the main part of the peak. If a giant pulse is recorded by detectors of 
low-energy charged particles which have low sensitivity to gamma-rays, there is no saturation prob- 
lem in most cases. However, severe difficulties are encountered instead when measuring the photon 
energy, obtaining spectra, and consequently deriving energy estimates. Nevertheless, the ent i re set 
of data on initial pulses, obtained both with gamma-ray spect rometers (jMazets et alJll979al . LL982J; 



Hur lev et al 



1999a; Mazcts et al. 1999a; Frederiks et al 



(Hurl ey et al.1 120051 : iPalmer et al.1 120051 : iTerasawa et al 



2005 



2007a|) and with charged-particle detectors 



Schwartz et al.l 120051 : iTanaka et al 



20071 ). can be used to establish the general pattern of the initial pulse. 



We will restrict this comparison of GFs to a limited number of their characteristics, for the 
following reasons. The time structure of the initial pulse is likel y to be very compl ex. The profile 
of both the steep leading edge of the pulse, as in SGR 18 06-20 (IPalmer et al.ll2005l ). and its more 
prolonged decay, as in SGR 1900+14 (jTanaka et al.1 120071 ). contain significant intensity variations 
on both short and long timescales. Among other things, the appearance of the light curve will be 
strongly affected by the spectral variability of the emission. At high count rates close to saturation 
level, dead-time and pileup corrections of light curves and energy spectra are difficult, if not impos- 
sible. Therefore, we include in Table [3] only two fairly reliable characteristics of the initial pulse, 
tn and E p ; tR, the rise time, is the time between the detection of hard gamma-radiation above 
the background to the peak of the initial pulse, E p is the photon energy at the maximum of the 
time-integrated EF(E) spectrum of the pulse. Table [3] also contains data on the peak luminosity 
L max and the energy release Q of the initial pulses, as well as the energy in the GF tails Qtaii- 
In keeping with tradition, we suggest the same nomenclature for extragalactic and galactic SGRs, 
namely SGR followed by the coordinates (a, 5) of the center of the error box. 



- 6 - 



Table [3] shows that estimates of both the luminosities and energy outputs in GFs are spread 
over a range of several orders of magnitude. On the other hand, the values of Qtaii are practically 
equal. Thus it is particularly significant that the afterglow of GRB 070201 is similar to those of 
the other SGRs. A smal l spread in SGR Qta.ii values implies similar magnetic field strengths, from 



confinement arguments ([Thompson Duncan! Il995l ). The only exception so far is the absence of 



a tail in GF 980618 from SGR 1627-41, which is probably caused by the relative weakness of the 
magnetic field or its orientation. On the whole, Table 3 clearly demonstrates that both the temporal 
and energetic characteristics of the event on 2007 February 1 match the general pattern of a giant 
flare very closely. Beyond a doubt, we can conclude that this event is a GF which originated in 
SGR 0044+42 in M31. 

On the Russian side this work was supported by Federal Space Agency of Russia and RFBR 
grant 06-02-16070. KH is grateful for IPN support under NASA grants NNG06GE69G (the IN- 
TEGRAL guest investigator program) and NNX07AR71G (MESSENGER Participating Scientist 
program). 
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Fig. 1. — 17-1130 keV Konus-Wind light curve of GRB 070201 (dead time corrected) at the highest 
available time resolution (2 ms) presented on a linear count rate scale (a), and also on a logarithmic 
scale (b) to enhance the low-intensity portions. 
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T-T (s) 

Fig. 2.— Konus-Wind light curve in three energy bands: Gl (17-70 keV), G2 (70-300 keV), and 
G3 (300-1130 keV), and the hardness ratios G2/G1 and G3/G2. The vertical dashed lines indicate 
the four successive 64-ms intervals where the energy spectra were measured. 
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Fig. 3. — Deconvolved photon spectra of the burst accumulated over the (T-To) intervals 0-64 ms 
(a), 64-128 ms (b), 128-256 ms (c), and the time-integrated spectrum, 0-256 ms (d). 
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Fig. 4. — Waiting mode count rates N in the energy bands Gl (a) and G2 (b) averaged over 
94.208 s intervals (excluding the burst). For Gl, the first interval after To exhibits a count excess 
AiV of 4.3(j above the average background level. The bottom panel (c) shows how this excess AN 
varies if the beginning of the 94.208-interval is shifted by an interval AT forwards or backwards 
from T step by step, with AT=2.944 s. 
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Fig. 5. — UV image of the M31 galaxy ( Thilker et al.l 20051 ) and the 3a IPN error box of 
GRB 070201. 
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Table 1. Summary of spectral fits 



Time interval 


a 


E p X 2 /dof 


(s) 




(keV) 



0-0.064 0.52lg;i| 360±^ 31.6/35 

0.064-0.128 0.56±g;^ 128±^ 11.3/15 

0.128-0.256 1.06+°;^ 123+^ 19.3/16 

0-0.256 0.98t&ii 296t| 39.7/40 



Table 2. IPN error box of GRB 070201 





RA(2000) 


Dec(2000) 


Center 


00 h 44 m 32 s 


+42°14 / 21" 


Vertices 






1 


00 h 46 m 18 s 


+41°56'42" 


2 


00 h 41 m 51 s 


+42°52'08" 


3 


00 h 42 m 47 s 


+42°31'41" 


4 


00 h 47 m 14 s 


+41°35'35" 



Table 3. Some characteristics of known giant flares 





SGR 0526-66 a 


SGR 1627-41 b 


SGR 1900+14 c 


SGR 1806-20 d 


SGR 0952+69 e 


SGR 0044+42 f 




(LMC) 








(M81) 


(M31) 


Data, YYMMDD 


790305 


980618 


980827 


041227 


051103 


070201 


Distance, kpc 


50 


10 


15 


15 


3600 


780 


Rise time, tji, ms 


~ 25 


~ 15 


~ 15 


~ 25 


~ 6 


~ 20 


Peak energy, E p , keV 


- 500 


- 150 


> 250 


- 850 


~ 900 


- 300 


Luminosity, -L ma x, erg s 


6.5 x 10 45 


3.4 x 10 44 


2.3 x 10 46 


3.5 x 10 47 


4.3 x 10 48 


1.2 x 10 47 


Energy release, Q, erg 


7 x 10 44 


1 x 10 43 


4.3 x 10 44 


2.3 x 10 46 


7 x 10 46 


1.5 x 10 45 


Tail energy, Q tail , erg 


3.6 x 10 44 


absent 


1.2 x 10 44 


2.1 x 10 44 


not detected 


7 x 10 43 



a The previous dead-time correction for the event l|Mazets et al.lll979al. fl982h has been revised here, assuming that the initial pulse 
displays the expected strong hard-to-soft spectral evolution. This resulted in a reduction in tn and an increase in L max and Q. 

iMazets etall l ll999ar ); IXptekar et all [bOOll l. 

kanaka etldl ^OOTj l ; iHurlev et al.l h999al l; IMazets et alj il999bll . 
^Palmer et"al] ll2005h ; iFrederiks et all j2007ah . 
HFrederiks et al.l j2007bh . 
f This work. 



